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Consider a reflecting diffusion in a domain in M.'^ that acquires drift in proportion to the 
amount of local time spent on the boundary of the domain. We show that the stationary 
distribution for the joint law of the position of the reflecting process and the value of the drift 
vector has a product form. Moreover, the first component is the symmetrizing measure on the 
domain for the reflecting diffusion without inert drift, and the second component has a Gaussian 
distribution. We also consider processes where the drift is given in terms of the gradient of a 
potential. 

1 Introduction 

This article is concerned with the higher dimensional version of a one-dimensional model originally 
introduced by Knight [21] and studied in more detail in [34, 35]. Computer simulations presented 
in [10] led to the conjecture that the stationary distribution for this higher dimensional process has 
a certain interesting structure. We prove this conjecture, and moreover answer questions about the 
stationary distribution left open in [10]. 

We start with a presentation of the model in a simple case, that of [10]. We consider a bounded 
smooth domain D C M.'^ and reflecting Brownian motion Xt in D with drift Kf. Let Bt be d- 
dimensional Brownian motion, n(x) be the unit inward normal vector of at x G dD and let Lt be 
the local time of X on dD, that is, a nondecreasing one-dimensional process with continuous paths 
that increases only when X^ G dD. The pair of processes {X, K) has the following representation: 



JO 

with {Xt,Kt) D X for all t > 0. Note that Xt is reflecting Brownian motion in D with 
normal reflection at the boundary and with drift Kt, Kt is an M'^-valued process that represents 
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the accumulated local time on the boundary in the direction normal to the boundary, and the drift 
Kf does not change when Xt is in the interior of D. We call the process K "inert drift" becaTise it 
plays a role analogous to the inert particle in Knight's original model [21] in one dimension. The 
main simulation result of [10] suggests that a stationary distribution for {X, K) exists and has a 
product form, i.e., that Xt and Kt are independent for each time t under the stationary distribution. 
Moreover, the first component of the stationary distribution is the uniform probability measure on 
D. We prove rigorously in this paper that this indeed holds, and we further show that the second 
component of the stationary distribution is Gaussian. 

The product form of the stationary distribution was initially a mystery to us, especially since 
the components X and K of the vector (X, K) are not Markov processes. There are models known 
in mathematical physics where the stationary distribution of a Markov process has a product form 
although each component of the Markov process is not a Markov process itself. Examples may 
be found in Chapter VIII of [23], in particular. Theorem 2.1 on page 380. As we will see at the 
beginning of Section 4, the product form of the stationary distribution in our model comes naturally 
from a computation with infinitesimal generators. 

The Gaussian nature of the stationary distribution for K is already known in the one-dimension- 
al case; see [34, 11]. An interesting phenomenon observed in this paper is that if the inward normal 
vector field n in the equation for K is replaced by Tn for some constant symmetric positive definite 
matrix T, {X, K) continues to have a product form for the stationary distribution, but this time 
the component K has a Gaussian distribution with covariance matrix T. 

The main goal of this paper is to address the existence and uniqueness of the stationary distri- 
bution of normally reflecting Brownian motion with inert drift and to give an explicit formula for 
the stationary distribution. We also consider a larger class of reflecting diffusions, including what is 
sometimes known as distorted reflecting Brownian motion — see Theorems 5.2 and 6.2 below. Dis- 
torted reflected Brownian motion is the reflecting diffusion with generator ^ V(/9V) for a suitable 
function p. 

We start by showing in Section 2 the weak existence and weak uniqueness of solutions to an 
SDE representing a large family of diffusions with reflection and inert drift. More specifically, let 
p be a function on D that is bounded between two positive constants and A[x) = {aij{x))nxn 
be a matrix-valued function on R'* that is symmetric, uniformly positive definite, and each aij is 
bounded and on D. The vector u(x) := ^^(a;)n(a;) is called the conormal vector at x G dD. 



Let cr(x) = {aij{x)) be the positive symmetric square root of A{x). For notational convenience, we 
use di to denote For e C2(M'^), let 




(1.1) 
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Let b(x) be the vector whose fe*'* component is 

1 

^^^^^ " 2^^^'^^^^''"'^^^''''' 

Let B be standard d-dimensional Brownian motion and v a bounded measurable vector field on 
dD. Consider the following diffusion process X taking values in D such that for all t > 0, 

dXt = a{Xt) dBt + h{Xt) dt + M{Xt) dLt + Kt dt, 

t i-^ Lt is continuous and non-decreasing with Lt = lg£)(Xs) dLg, (1-2) 
dKt = viXt)dLt. 

In Theorem 2.1, we show that the above stochastic differential equation (SDE) has a unique weak 
solution {X,K) for every starting point (xo,fco) € D x R.'^. The solution {X,K) of (1.2) is called 
a (symmetric) reflecting diffusion on D with inert drift. Let X^ be symmetric reflecting diffusion 
on D with infinitesimal generator C in (1.1); that is, X^ is a continuous process taking values in D 
such that for every t>0, 



dXf = a{X^) dBt + h{X^) dt + u{X^) dL^t 



t H- > is continuous and non-decreasing with L° = IgniXg) dL^. 



(1.3) 



The continuous non-decreasing process L° is called the boundary local time of When a is the 
identity matrix, X° is distorted reflecting Brownian motion on D. The main observation of Section 
2 is that the reflecting diffusion {X, K) with inert drift can be obtained from the reflecting diffusion 
(^X^, Kq + Jp n{Xg ) dL^^ without inert drift by a suitable Girsanov transform, and vice versa. 

The questions of strong existence and strong uniqueness for solutions to (1.2) are discussed in 
Section 3. They are resolved positively under the additional assumption that v = oqU for some 
constant ao G M. This section uses some ideas and results from [24], but the main idea of our 
argument is different from the one in that paper, and we believe ours is somewhat simpler. 

In Section 4 we consider symmetric diffusions with drift given as the gradient of a potential. We 
do this because the analysis of the stationary distribution is much easier in the case of a smooth 
potential than the "singular" potential representing reflection on the boundary of a domain. More 
specifically, let F be a symmetric positive definite constant d x d-matrix and V G C^{D) tending 
to infinity in a suitable way as x approaches the boundary dD. Consider the following diffusion 
process X on D associated with generator C = ^e^ (e~^^V) but with an additional "inert" drift 
Kt: 

= aiXt) dBt + h{Xt) dt - i(AVF)(Xt) dt + Ktdt, 
= -lTVViXt)dt, 

Here A = A{x) = a^a is a d x d matrix-valued function that is uniformly elliptic and bounded, 
with aij G C^(D), and b = (6i,--- ,6^) with bk{x) := ^e^W (e-^(^)ajfc(x)) . We show 
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in Theorem 4.3 that if e G Wq''^{D), then the SDE (1.4) has a unique conservative solution 
{X, K) which has a stationary probabihty distribution 

Tr{dx,dy) = cilD(x)e-^(^) e'^^'^^^y^ dxdy. 

In other words, {X,K) has a product form stationary distribution, with ce~^ dx in the X com- 
ponent, and a Gaussian distribution with covariance matrix T in the K component. Observe 
that c dx with normahzing constant c > is the stationary distribution for the conservative 
symmetric diffusion X^ in D with generator ^e^V(Ve~^AV). The uniqueness of the stationary 
measure for solutions of (1.4) is addressed in Proposition 4.8 under much stronger conditions. 

In Section 5, we prove weak convergence of a sequence of solutions of (1.4) to symmetric 
reflecting diffusions with inert drift given by (1.2) where v = Fn for some symmetric positive 
definite constant matrix. This implies 

7r{dx,dy) = C21d{x)p{x) e"^^ ^'^''^^dxdy 

is a stationary distribution for solutions {X, K) of (1.2) with v = Fn; see Theorem 5.2. Observe that 
C3 p{x) dx, where C3 > is a normalizing constant, is the stationary distribution for the symmetric 
reflecting diffusion X^ on D with generator ^ V(/9^V) in (1.3). 

Finally, Section 6 completes our program by showing irreducibility in the sense of Harris for 
reflecting diffusions with inert drift given by (1.2) under the assumption that v = Fn for some 
symmetric positive definite constant matrix, A is the identity matrix, u = n, and b = ^Vlogp. 
Uniqueness of the stationary distribution follows from the irreducibility. 

The level of generality of our results varies throughout the paper, for technical reasons. We 
leave it as an open problem to prove a statement analogous to Theorem 6.2 in the general setting of 
Theorem 2.1. The calculation at the beginning of Section 4 indicates that in order for the solution 
(X, K) of (1.2) to have a product form stationary distribution, the inert drift vector field v has to 
be of the form F n for some symmetric positive definite constant matrix F. We leave its verification 
as another open problem. 

Our model belongs to a family of processes with "reinforcement" surveyed by Pcmantlc in [30]; 
see especially Section 6 of that survey and references therein. Papers [5, 6, 7] study a process with 
a drift defined in terms of a "potential" and the "normalized" occupation measure. While there is 
no direct relationship to our results, there are clear similarities between that model and ours. 

We are grateful to Boris Rozovsky and David White for very useful advice. 

2 Weak existence and uniqueness 

This section is devoted to weak existence and uniqueness of solutions to an SDE representing a 
family of reflecting diffusions with inert drift. 
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Let D he a bounded domain in M'^, d > 2, and denote by n the inward unit normal vector 
field on dD. Throughout this paper, all vectors are column vectors. Let p be a function on 
D that is bounded between two positive constants and A{x) = {aij{x))nxn be a matrix-valued 
function on that is symmetric, uniformly positive definite, and each is bounded and on 
D. The vector u(x) := ^A{x)n{x) is called the conormal vector at x e dD. Clearly there exists 
ci > such that u{x) ■ n(x) > ci for all x G dD. Let a{x) = {aij{x)) be the positive symmetric 
square root of A{x). For (p G C^{W^), let 



1 

2^ {p{x)aij{x)djip{x)) . 



Let h{x) be the vector whose A;*'* component is 
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Thus is the same as C operating on the function fkix) = Xk- 

Let B be standard d-dimensional Brownian motion and v a bounded measurable vector field 
on dD. Consider the following system of stochastic differential equations, with the extra condition 
that XtGD for all t > 0: 

dXt = a{Xt) dBt + h{Xt) dt + u{Xt) dLt + Kt dt, 

1 1— > Li is continuous and non-decreasing with Lt = Jg 1qd{Xs) dLg, (2-1) 
dKt = Y{Xt)dLt. 

The proof of the next theorem says that reflecting diffusions with inert drift can be obtained 
from the corresponding symmetric reflecting diffusions without inert drift by suitable Girsanov 
transforms, and vice versa. 

Theorem 2.1 For every x £ D and y G M*^ there exists a unique weak solution {{Xt, Kt),t G 
[0,oo)} to (2.1) with {Xq,Kq) = {x,y). 

Proof. Consider the following SDE, 

(dXt = a{Xt) dBt + hiXt) dt + u{Xt) dLt, 



t. ...... (2-2) 



1 1 1— >■ is continuous and non-decreasing with Lt = Jq ldD{Xs) dLg, 

with Xt E D for every t > 0. Weak existence and uniqueness of solutions to (2.2) follow from [12] 
or [15]. (We have been informed by the authors of [15] that there is a gap in the proof of Case 2 in 
[15], but we need only Case 1. Moreover they have informed us that a correction for Case 2 is in 
press.) The distribution of the solution to (2.2) with Xq = x & D will be denoted by P^;. 
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Note that the remaining part of the proof uses only the C -smoothness of the domain and 
Lipschitz continuity of aij and p. Let Kt := y + Jq v{Xs) dLg and a~^{x) be the inverse matrix of 
a{x). Define for t > 0, 

Mt = exp ^* (T-\Xs)Ks dBs -^j^ \a-\Xs)Ksf ds^ . 

It is clear that M is a continuous positive local martingale with respect to the minimal augmented 
filtration {Tt,t > 0} of X. Let T„ = inf{t > : \Kt\ > 2"}. Since Lt < oo for every t < oo, 
Pj^-a.s., and jvj is uniformly bounded, we see that Kf < oo for every t < oo, Pa;-a.s. Hence, 
Too '■= liiUn^oo^n = CO) P^-a-s. For every n > 1, {MTn/^tjJ^tj't > 0} is a martingale. 
For X e D, y eW^ and n > 1, define a new probability measure Qx,y by 

dQx,y = Mt„ dFx on J^t^ for every n > 1. 

It is routine to check this defines a probability measure Qx,y on J^^o- By the Girsanov theorem (cf. 
[31]), the process 

Wt:=Bt- I a-\X,)K,ds, 
Jo 

is a Brownian motion up to time T„ for every n > 1, under the measure Qx,j/- Thus we have from 
(2.2) that under Qx,y, up to time for every n > 1, 

dXt = a{Xt) dWt + h{Xt) dt + u{Xt) dLt + Kt dt. 

In other words, {{Xt, Kt),0 <t< T^o} under the measure Qx,y is a weak solution of (2.1). 

We make a digression on the use of the strong Markov property. At this point in the proof, we 
cannot claim that {X, K) is a strong Markov process. Note however that if T is a finite stopping 
time, then [Xt+T ■, Ki+t) is again a solution of (2.1) with initial values {Xt,Kt)- We can therefore 
use regular conditional probabilities as a technical substitute for the strong Markov property; this 
technique has been described in great detail in Remark 2.1 of [3], where we used the name "pseudo- 
strong Markov property." Throughout the remainder of this proof we will use the pseudo-strong 
Markov property in place of the traditional strong Markov property and refer the reader to [3] for 
details. 

We will next show that Too = oo, Q^^.jz-a.s., i.e., the process is conservative. This is the same 
as saying \K\ does not "explode" in finite time under Qx,y The intuitive reason why this should 
be true is the following. Consider a one-dimensional Brownian motion starting at 1 with a very 
large constant negative drift of size c and reflect it at the origin. Then a simple calculation shows 
that the local time accumulated at the origin by this process up to time 1 is also of order c. This 
suggests that if K is very large, the local time accumulated in a time interval of order 1 by X on 
the boundary dD will be approximately proportional to K. Since this feeds back into the right 
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hand side of the definition of K in (2.1), one would expect K to grow at most exponentially fast 
in time. 

Consider £ = 2^^ > where j > 1 is an integer. Our argument applies only to small e > so 
we will now impose some assumptions on e. Consider xq G dD and let CSxq be an orthonormal 
coordinate system such that xq = in CSxq and the positive part of the d-th axis contains n(xo). Let 
no = n{xo). Recall that D has a boundary and that there exists ci > such that u(x)-n(x) > ci 
for all X G dD. Hence there exist £o > ^'^d C2 > 0, such that for all e G (0, Sq), every xq G dD 
and all points x = {xi, . . . ,X(i) G dD n B{xq, (6/c2 + 5)e), we have \X(i\ < e/2 and u(x) ■ no > C2, 
in CSxQ- Since |v(x)| < C3 < oo for all x G dD, we can make £0 > smaller, if necessary, so that 
(2c2)/(c3£) - 5e > e for all e G (0,eo]- 

Let 

So = M{t > : \Kt\ > 1/e}, 
Sn+i = mf{t >Sn:\ \Kt\ - \KsJ \ > {603/02)8}, n>0. 

We will estimate Qx,y{Sn+i - Sn > £^ \ J^s,,) for n = 0, . . . , [e"2c2/(6c3)]. Note that Xs„ G dD 
for every n because K does not change when X is in the interior of the domain. Note also that 
\Ks„ \ < 2/e for every n < £"^02/(603) by construction. 
For n > 0, let 



, , rSn+t rSn+t 

/ a{Xs)dWs+ / h{Xs)ds, 



Fn=\ sup <el. 

[te[0,e^] J 

It is standard to show that there exist £o,Po > 0, not depending on n, such that if £ < £0, then 

Qx,y{Fn\J'sJ>PO- 

Let 

Rn = {Sn + £') A inf{t > Sn : \Kt\ > 4/e}. 

Suppose that the event Fn holds. We will analyze the path of the process {{Xt, Kt), Sn < t < 
Rn}- First, we will argue that Xt G B{Xs„, (6/c2 + 5)£) for Sn <t < Rn- Suppose otherwise. Let 
Ui = inf{i > Sn-.Xti B{Xs^, (6/c2 + 5)£)} and U2 = sup{t < Ui : Xt e dD}. We have 



{6/02 + 5)e = \Xu,-XsJ 
So on Fn n {Ui < Rn}, 



rU2 rUi 

/ u{Xt)dLt = / u{X, 

" Sfi J Sn. 



Yu- 



Ui 



,-s^+ Ktdt+ u{Xt) 

*^ Sn *^ Sn 



dLi 



t)dLt 
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> (6/c2 + 5)£ 



rUi 

l-Sn + / 



Ktdt 



> (6/c2 + 5)£ - £ - / 

Js 

> (6/c2 + 5)£ - £ - 4£ 
= (6/c2)£. 



4/edt 



We will use the coordinate system CSxg^ to make the following observations. The last formula 
implies that the d-th. coordinate of Jg^ ^{^t) dL^ is not less than 6e. Hence the d-th coordinate of 
Xjj^ must be greater than or equal to 



6s- 



Yui-Sr^ + Ktdt >6£-e- A/edt > 6e - e - 4e = e. 

" Sn. Sn. 



This is a contradiction because G dD and the d-th coordinate for all x G dD f] B{Xs„, (6/c2 + 
5)£) is bounded by e/2. So on Fn, we have Xf G B{Xs„, (6/c2 + 5)e) for t € [Sn, Rn]- 

We will use a similar argument to show that Rn = Sn + £^ on Fn- Suppose that Rn < Sn + ■ 
Then Sn < Rn and \Kji^\ = 4/e. Let Us = sup{t < Rn '■ Xt E dD}. The definition of Sn implies 
that Ks„ < 2/£ forn < [£-202/(603)]. We have 



/ 

JSr 



Y{Xt)dLt 



\KR^-KsJ>4/s-2/s = 2/e. 



Since |v(a;)| < C3, we have Lr^ — Lg^ > 2/(c3e), so the d-th coordinate of fg " u{Xt) dLt, which is 
the same as Jg^ u(Xt) dLt, is bounded below by (2c2)/(c3e). But then on F„, the d-th coordinate 
of Xug must be greater than or equal to 



(2C2)/(C3£) - 



YRr. 



rRn rSn+e^ 
-5„+ / Ktdt > (2C2)/(C3£) - £ - / 4/£dt= (2c2)/(c3£)-5£>£. 

■J Sn. 'J Sn 



This is a contradiction because Xu^ G dD and the d-th coordinate for all x G dDCiB^xo, (6/c2-|-5)£) 
is bounded by e/2. So on Fn we have i?„ = Sn + 

We will now use the same idea to show that on F„, -Z^^^+^a — Ls^ < (6/c2)£. Assume that 
^Sn+e'^ — Ls„ > (6/c2)£. Let U4 = swp{t < Sn + £^ '■ Xt E dD}. The d-th coordinate of 
Jg'^ u{Xt) dLt is bounded below by 6£. But then on F„, the d-th coordinate of Xu^ must be greater 
than or equal to 



\yRn-Sn + 



Ktdt 



Sn 



> 6e 



4/edt > e. 



This is a contradiction because Xu^ G dD and the d-th coordinate for all x G dDnB{xo, (6/c2-|-5)£) 
is bounded by £/2. We see that if the event Fn holds, then 



\K. 



-\Ksj\< r 

JSr 



\v{Xt)\ dLt < C3{Ls„+e^ - LsJ < C3(6/C2)£. 
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Hence, if the event F„ holds, then Sn+i > + e^. We see that 

Recall that wc took e = 2^^ so that = Tj. The last estimate and the pseudo-strong Markov 
property applied at stopping times Sn allow us to apply some estimates known for a Bernoulli 
sequence with success probability po to the sequence of events {Sn+i > S'n + £^}- Specifically, there 
is some pi > so that for all sufficiently small e = > 0, 

Qx,y{Tj+l - Tj > (c2/(6c3))po/2 I J^T,) 

'[£-2c2/{6c3)] 



^ {Sn+1 - Sn) > (c2/(6c3))po/2 



n=0 



2 [£ 2C2/(6C3)] 



> Pi- 



n=0 



Once again, we use an argument based on comparison with a Bernoulli sequence, this time with 
success probability pi. We conclude that there are infinitely many n such that Tj+i — Tj > 
(c2/(6c3))po/2, Qx,y-^-S- We conclude that Too = oo, Qa;,j,-a.s., so our process {Xt,Kt) is defined 
for ah t G [0, oo). 

Next we will prove weak uniqueness. Suppose that Q'x^y is the distribution of any weak solution 
to (2.1) and define by 



^x = j^ dQx,y on JP-T„ for every n > 1. 



Reversing the argument in the first part of the proof, we conclude that X under solves (2.2). 
It follows from the strong uniqueness for (2.2) that = Px, and, therefore, Q^,?/ ~ Qx,y on J-Tn- 
Since this holds for all n, we see that Q'x^y = Qx,y on J^^. □ 



Remark 2.2 The assumptions that D is a bounded domain and that the Ojj's and p arc on 
D are only used in the weak existence and uniqueness of the solution X to (2.2). The remaining 
proof only requires that D be a bounded domain and that the aij's and p are Lipschitz in D. In 
fact, when D is a bounded domain and the 's and p are Lipschitz on D, the weak existence of 
solutions to (2.2) follows from [12] and so we have the weak existence to the SDE (2.1). We believe 
the weak uniqueness for solutions to (2.2) and consequently to (2.1) also holds under this weaker 
assumption by an argument analogous to that in [2, Section 4]. However to give the full details of 
the proof would take a significant number of pages, so we leave the details to the reader. □ 
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We remarked earlier that if T is a finite stopping time, then (Xt+T, Kt-^-r) is again a sohition 
to (2.1) with starting point {Xt,Kt)- This observation togctlicr witli the weak uniqueness of the 
solution to (2.1) implies that (X^jK^) is a strong Markov process in the usual sense; of. [1, Section 



3 Pathwise uniqueness 

We will prove strong existence and strong uniqueness for solutions to (2.1) under assumptions 
stronger than those in Section 2, namely, we will assume that the vector field v is a fixed constant 
multiple of u. Throughout this section, D is a bounded domain in M'^, each aij and p are C^'^ 
on D and so the vector b in (2.1) is Lipschitz continuous. Our approach to the strong existence 
and uniqueness for solutions to (2.1) uses some ideas and results from [24], but the main idea of 
our argument is different from the one in that paper, and we believe ours is somewhat simpler. It is 
probably possible to produce a proof along the lines of [24] , but a detailed version of that argument 
adapted to our setting would be at least as long as the one we give here. 

Theorem 3.1 Suppose that v = qqu for some constant oq G M. For each {x,y) G Z) x R"', there 
exists a unique strong solution {{Xt, Kt),t G [0,oo)} to (2.1) with {Xq,Kq) = {x,y). 

Proof. When = 0, Kt = Kq for every t > 0. In this case, the result follows from [15] (as 
mentioned above there is a gap in the proof of Case 2 in [15]; we need only Case 1). So without 
loss of generality, we assume ao ^ 0. 

First we will prove pathwise uniqueness. Suppose that there exist two solutions (Xt,Kt) and 
(X^, K^), driven by the same Brownian motion B, starting with the same initial values {Xq, Kq) = 
{Xq, Kq) = {x,y), and such that {Xt, Kt) / {X^, K[) for some t with positive probability. 

According to Lemma 4.1 of [24], there exists a matrix-valued function A(x) = {Kj{x)}\<ij<d, 
X G R'^, such that x — A(x) is uniformly elliptic and in C| and such that the following formulas 



1.5]. 



(3.1) and (3.2) hold. We have 



u(x)^A(x) = n{x)^ for X G dD. 



(3.1) 



Moreover, there exists ci < oo such that 



C]\x - x'\^ + u(x)'^A(x)(x' - x)>Q for X e dD and x' G D. 



(3.2) 



In our setting, we can take A = 2^ . Let C2 > 1 be such that 



C2 Idxd < A(x) < C2ldxd 



for every x G M' 



d 



(3.3) 



where Idxd is the d x d-dimensional identity matrix. Define A := cf. 
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Let 



Ut = f a{Xt) dBt + /* h{Xt) dt, 
Jo Jo 



Vt= f Ktdt, 
Jo 

and define U[ and V{ in a similar way relative to X' . 

We fix an arbitrary pi < 1 and an integer fco such that in fco Bernoulli trials with success 
probability 1/2, at least A;o/4 of them will occur with probability p\ or greater. Let t\ = 1/ko. 

Consider some e, cq > 0, and let 

Ti = inf{t > O:\Xt- > or \Ut - Ul\ > or \Kt - K[\ > 0}, 
n = {Tk-1 + h) A inf{t > Tk-i : \Xt - X^\ V \Ut - Ui\ > X'^-'e} 

A inf{i > Tk-i : \Lt - Lt,_, | V |L; - L^^_^ | > cq}, k > 2. 

We will specify the values of e and cq later in the proof. 

It is easy to see from (2.1) that it is impossible to have Xt equal to X^ on some random interval 
[0,T*], and at the same time Ut / C// or Kt / for some t G [0,T*]. Hence, with probability 1, 
for every 5 > there exists G (Ti ,Ti -\- S] such that Xt^ ^ X[ . 

The idea of our pathwise uniqueness proof is as follows. Define 

5 = (Ti + 1/4) A inf{t > Ti : (L* - LtJ V (L^ - L^J > co/(4ti)}. 
We will show that for any pi <1 and integer Uq = ko{pi) as defined above and every £ > 0, 

P ( sup \Xt - X'tl < eX'"' ] > pi. 
\te[o,S] J 

As e > is arbitrary, we have 

P {Xt = X[ for every te[Q,S])>pi> 0, 

which contradicts the definition of Ti , in view of the remarks following the definition. This contra- 
diction proves the pathwise uniqueness. 

Gronwall's inequality says that if g{t) is nonnegative and locally bounded and g{t) < a + 
b jQg{s)ds, then g{t) < ae^*. Suppose now that / is a nonnegative nondccrcasing function and 
g{t) < f{t) + b Jq g{s) ds. Applying Gronwall's inequality for t <ti with a = f{ti), we have 

9{ti)<f{h)e^K (3.4) 

We apply the inequality with g{t) = \Kt — K[\ and 

f{t) = \KT,-K'T^\+ sup [\Xs-X',\ + \Us-U',\). 

Ti<s<t 
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Since v = oqU, we have Jq u{Xs) dLg = -^{Kt — Kq) and so 

Kt-Ko = aoiXt -xo-Ut- Vt), 

and similarly for K'. Thus 

Kt - K[ = ao(Xt - X[ -Ut + U[ -Vt + F/), 

and hence for t >Ti, 

\Kt - K[\ < \Kt, -K!r,\ + 2\ao\ I sup - X'J + |C/, - U'J) + f \Ks - K'^\ ds] . 

\Ti<s<t JTi J 

By (3.4), fort>ri, 

\Kt - K[\ < e2|«ol(*-^i) (\Kt, -K^,\+ 2\ao\ sup - X',] + |C/, - C/^|) ) . 

Y ri<s<t J 

Recall that h = l/ko- It follows that - Ti < 1 and e^l^oK*-^!) < e^l^d for t < T^^. The 
definition of the T^'s implies that 

sup {\Xs - X'^\ + 1^7, - U',\) < 2X''-h. 

0<s<Tk 

We obtain for t <Tk with k < ko, 

\Kt - Ki\ < e2|«ol(l + 4|ao|)A^-^£. (3.5) 

We have 

S := {Xt, - X'^f K{Xt,){Xt, - J - {Xt,_, - X'^^_f K{Xt,_,){Xt,_, - 
= r d{{Xt-XifA{Xt){Xt-Xi)) 

<C3 f ' A{Xt){Xt - Xi) ■ {{a{Xt) - a{Xi)) dBt 

+ C4 / ' K{Xt){Xt-X[)-{h{Xt)-h{X[))dt 

+ C5 f ' A{Xt){Xt - Xi) ■ {Kt - K[) dt 

+ f ' K{Xt){Xt-X[)-n{Xt)dLt 

- c, r K{Xt){Xt-X[).n{X't)dL't 
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+ Cg 



/ ' \Xt- dBt + h{Xt) dt + Ktdt + vi{Xt) dLt) 

JTk-i 



Rewrite the second to the last term as 

(■Tk 



-C7 r A{X^)iXt-X',).n{X[)dL',-C7 f" {K{Xt) - KX't)){Xt - X'^ ■ n{X'^ dL[. 



We apply (3.5) to J^^_^ A{Xt)iXt - X'^.) ■ {Kt - K[) dt and use the Lipschitz property of A to 
obtain 



E<cs [ ' A{Xt){Xt - Xi) ■ {a{Xt) - a{Xi)) dBt 

+ cj ' A{Xt){Xt - X[) ■ {h{Xt) - h{Xi)) dt 
+ C9 r \K{Xt){Xt-Xi)\\^-^edt 



+ CIO / \Xt-Xl\ULt 
iTk-i 



+ cu f ' \Xt-XifdL[+ ci2 / ' \Xt-XlfdL[ 

[ ' \Xt- XiWaiXt) dBt + hiXt) dt + Ktdt + u(Xt) dLt) 



+ Cg 



(3.6) 



Here, we made use of (3.2) to bound the two terms involving integrations with respect to the local 
times Lt and L[. 

It follows from (3.3) that there exists C13 > such that if 

(XT,-X^J^A(XTj(XT,-X^J-(Xr,_,-X^^_J^A(Xr,_J(XT,_,-X^,_J < ciaA^Me^ (3.7) 

and 



then 



In view of our assumptions on A and cr, 



(3.8) 
(3.9) 



Var (c3 / ' A{Xt){Xt - X't) ■ {a{Xt) - a{X^)) dBt I J't.A 
< E (ci4 / ' sup \Xs - X'f dt Jr,_i) 
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Wc make ti > smaller (and therefore ko larger), if necessary, so that by Doob's and Chebyshev's 
inequalities, 

C3 r A{Xt){Xt-X'^)-{a{Xt)-a{X[))dBt > (1/10)ci3A2M£2 J < 1/8. (3.10) 



We make > smaller, if necessary, so that 



C4 f' K{Xt){Xt - X't) ■ {h{Xt) - b(X;)) dt < ci5 r 



< C15 / \Xt-Xi\^dt 
1 



(3.11) 



/ sup \X,-X'fdt 



< Ci5 

< Ci5tlA2('=-l)£2 < (1/10)ci3A2('=-1)£^ 



If ti > is sufficiently small then 



C9 / |A(XO(Xi-X;)|A^"i£dt 



< ci7 / sup -X^IA^'^^edt (3.12) 



< (1/10)ci3A2(^-i)£^ 



We make cq > in the definition of so small that 

< (l/10)<;,3A2<'-'>e^ 



r 



cio |Xt-X;pdLt + cii /"^' \Xt-XifdL[ + ci2 f^' \Xt-X[\^dL[ 



(3.13) 



Completely analogous estimates show that for sufficiently small ti,co > 0, 



C8 / \Xt - X[\\a{Xt) dBt + h{Xt) dt + Ktdt + u{Xt) dLt] 



■ r 



> (1/10)ci3A2('=-i)£2 < 1/8. 



(3.14) 



Combining (3.6) with (3.10)-(3.14), we see that conditioned on J^T^-n with probability 3/4 or 
greater the following event holds: 

{Xt, - X'TfK{XT,){XT, - X't^) - {Xt,_, - X'T,_fk{XT,_,){XT,_, - J < ci3A2(^-i)£2. 

In view of (3.7)-(3.9), this implies that if (3.8) holds, then 

\XT,-X'T,^\<\^-^e. (3.15) 

Let 

Tl = {Tk-i + h) A inf{t > Tfe_i : \Xt - Xi\ > X'^'^e}. 
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By Doob's inequality, if ii > is sufficiently small, 

F ( sup \Us -U',\> (l/2)A*^-ie J'n A < 1/4. (3.16) 

\Tk-i<s<T;^ J 

Let 

Fk = {Tk = Tk-i + h} U {\Lt, - Lt,_,\ V \L't^ - L't^J > cq}. 
By (3.15) and (3.16), 

n^fc|^T,_J>l/2. (3.17) 

Recall the definition of relative to pi G (0, 1). Repeated application of the strong Markov 
property at the stopping times Tk and comparison with a Bernoulli sequence prove that at least 
/co/4 of the events Fk will occur with probability pi or greater. This implies that with probability 
pi or greater, 

either T^, - Ti > 1/4 or Lt,^ - Lt, > co/(4ti) or L^^^ - L'^^ > co/{Ui). 
Prom the definitions of T\ and 5, 

P f sup \Xt - X[\ < eX^° ) > pi. 
V*G[o,s] ) 

As e > is arbitrary, we have 

P {Xt = X[ for every i G [0, 5]) > pi > 0, 

which contradicts the definition of Ti. This completes the proof of pathwise uniqueness. 

Strong existence follows from weak existence and pathwise uniqueness using a standard argu- 
ment; see [31, Section IX. 1]. n 



4 Diffusions with gradient drifts 

This section is devoted to analysis of a diffusion with inert drift given as the gradient of a potential. 
We will use such diffusions to approximate diffusions with reflection. The analysis of the stationary 
measure is easier in the case when the drift is smooth. Let F be a function on D that goes 
to +00 sufficiently fast as x approaches the boundary of D. We will consider the diffusion process 
on D associated with generator C = V {e~^ AV) but with an additional "inert" drift. More 
precisely, let T be a non-degenerate constant d x d-matrix. We consider the SDE 

dXt = aiXt) dBt + h{Xt) dt-\ {AW) {Xt) dt + Ktdt, 

dKt =-lrvv{Xt)dt, 
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where A = A{x) = {aij) = a^a is uniformly elliptic and bounded. We assume that a € C^{D) 
and b = {hi,--- ,bd) with bkix) := ^e^(^) X^li (e"'^^''^aifc(a;)) • Note that, since a G C^(D), 
V G C^{D) and b is bounded, (4.1) has a unique weak solution {X,K) up to the time inf{t > : 
Xt^D 01 Kt = oo}. 

To find a candidate for the stationary distribution for {X, K) , we will do some computations 
with processes of the form f{Xt,Kt), where / G C'^{M?'^). We will use the following notation, 

Vj;/(x,y) 
^xf{x,y) 

For / G C2(M2d)^ by Ito's formula, we have 

1 

df{Xt, Kt) = VxfdXt + Vyf dKt + - ^ didjfd{X\ X^)t 

= local martingale + (^xf + ^xf ■ Kt — ^^yf " rVj;^^ dt. 
So the process {X, K) has the generator 

y) := £xf{x, y)+y V,/(ar, y) - V.F (x) • Vyf{x, y) (4.2) 

for a; G -D and y G M*^. 

We will now assume that (X, JT) has a stationary measure of a special form. Then we will do 
some calculations to find an explicit formula for the stationary measure, and finally we will show 
that the calculations can be traced back to complete the proof that the measure we started with is 
indeed the stationary distribution. 

Suppose that [X^K) has a stationary distribution vr of the form Tr{dx,dy) = pi{x) p2{y) dx dy . 
It follows that Q*t: = 0, in the sense that for every / G Dom(^), 

/ Qf{x, y)Tr{dx, dy) = 0. 
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d d 
— /(xi , . . . , 2/1 , . . . , yd) , . . . , -Q^Jixi ,...,Xd,yi,...,yd)] , 



T 



d d 

—f{xi, . . . ,Xd,yi, . . . ,yd), ■ ■ ■ , -g;^fixi, . . . ,Xd,yi, . . . ,yd) ) , 



dxi 



iE|-(e-''«««w4(e''«/(x,.)). 



Then we have for every f E C^{D x M"), by the integration by parts formula, 

= J^^(^J^ipiix)y ■Vxf{x,y) + pi{x)£xf{x,y)) dx^ P2{y)dy 
~lJj,{L, P2(y)r V,.F(a;) • Vyfix, y) dy^ pi(x) dx 
= 1(1 {-^xPi{x)-yf{x,y) + Clpi{x)f{x,y)) dx] p2{y)dy 
\ (^J^^ r V^Vix) ■ VyP2{y)f{x, y) dy^ pi{x) dx . 



+ 2 



This implies that 



P2{y) (V^piix) ■ y - Clpi{x)) - ]^pi{x)TVxV{x) ■ VyP2{y) = for e x R<^. (4.3) 

We now make an extra assumption that pi (x) = ce~^(^) . Then we have 

P2{y)VxV{x) ■y + \T VxV{x) ■ VyP2{y) = for {x, y)eDx R'^. 

Since V{x) blows up as x approaches the boundary dD, {VxV{x),x G D} spans the whole of R'^. 
(If not, there exists v eBJ^ such that {v, W{x)) < for every x G M*^. One then gets a contradiction 
to the fact that there exists ro > (possibly ro = +oo) such that limr-^-r-o V{x + vr) = oo for every 
X G D.) So we must have 

y + log p2 (y) = for every y G M'*. 

This implies that F is symmetric and that 

Vj, log p2{y) = -2T~^y for every y G R'^. 

Hence we have 

logP2(y) = -{T~^y,y) + ci, 



or 



P2(y) = C2exp ( - (r ^y,y)y 



The above calculations suggests that when F is a symmetric positive definite matrix, the sta- 
tionary distribution for the process {X,K) in (2.1) has the form 



C3 1d{x) exp ^ - V{x) - (F ^y, yj^ dx dy, 



where C3 > is the normalizing constant. This is made rigorous in the following result. Recall that 
a process is said to be conservative if it does not explode in finite time. In the present setting, since 
we consider processes taking values on D x W^, the is means that the second component does not 
explode and that the first component does not reach the boundary of D. 
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Theorem 4.1 Suppose that T is a symmetric positive definite matrix, each CTij is C on D so that 
A = a^a is uniformly elliptic and bounded on D, and V is a potential on D. Suppose that 

Tr{dx,dy) := co1d{x) exp ^ - V{x) - {r~^y,yfj dxdy 

is a probability measure on D xW^ such that the dijfusion process of (4.1) with initial distribution vr 
is conservative. Then the process of (4.1) has ir as a (possibly not unique) stationary distribution. 

Proof. Let g be the operator defined by (4.2) with domain V{g) = C^{D x W^). The above 
calculation shows that 

[ gf{x,y)Tridx,dy) =0 for every / G P(a). (4.4) 

Let E := D xM.'^ and = ELl{A} be the one-point compactification of E. Recall g is defined 
by (4.2) with V{g) = Cl{E). Since (4.4) holds, we have by Theorem 9.17 in Chapter 4 of [16] 
that TT is a stationary measure for some solution P to the martingale problem for {g,T){g)). In 
[16, Theorem 4.9.17], the measure P is a probability measure on the product space However 
by [16, Corollary 4.3.7], any solution of the martingale problem for (^,P(^)) has a modification 
with sample paths in the Skorokhod space D(M+, E^) of right continuous paths on E a having left 
limits. Thus we can assume that P is a ©(M^., £^A)-solution to the martingale problem {g,T>{g)) 
with stationary distribution tt. Let (X, K) denote the coordinate maps on E^ and set 

C:=inf{t >0: (Xt,i^O = A}. 

We show next that {{Xt, Kt), t < C} under P is a solution to SDE (4.1) up to time C,. 
Let f{x,y) = g{x)h{y) with g G Cl{D) and h G Cl{^'^). Then under P, 

g{Xt)h{Kt) = giXo)hiKo) + martingale (4.5) 

+ ^ (h{Ks)CJ{Xs) + h{Ks)Ks ■ V:cg{Xs) 

- ^g{Xs)VV^V{Xs) ■ Vyh{Ks)) ds. 

Here and below, we use the convention that for any function / on E, /(A) := 0. Since tt is the 
stationary distribution for [X,K), by the dominated convergence theorem, the above holds for 
g G C'^{D) and h G C'^{W^) that have bounded first derivatives. So in particular, we have for 
g G C^iD), under P 

g{Xt) = g{Xo) + martingale + / [jOxgiXs) + Ks ■ Va,g{Xs)) ds for t > 0. 

Jo 

Let g{x) = Xi. There are local martingales M = (M^, • • • , M'^) such that 

dXt = dMt + h{Xt) dt-^ (AW) {Xt) dt + Kt dt, t < (. 
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Applying Ito's formula to XiXf yields {M\AP)t = {X\X^)t = j*^ aij{Xs) ds for t < (. Define 
dBt := a^^{Xt) dMf. Then {Bt, t < (} is a. Brownian motion up to time ^ and Mt = Jq (t{Xs) dBg 
lor t < Hence we have 

dXt = a{Xt)dBt + h{Xt)dt-^{AVV){Xt)dt + Ktdt, t<C. (4.6) 

Prom (4.5), we have for every h G C^(R'') that has bounded derivatives, 

1 /•* 

h{Kt) = h{Ko) + local martingale - - / TVxV{Xs) ■ Vyh{Ks) ds, t <(. 
In particular, taking h{y) = yi, I < i < d, there are local martingales A*" = {N^, • • • , N'^) such that 

dKt = dNt - ^rVxV{X,) ds, t < c. 
Applying Ito's formula to K^K^ yields {N\ N^)t = {K\ K^)t = for t < C- Hence 

dKt = -]^VVxV{Xs)ds, t<C. 

This together with (4.4) implies that {{Xt,Kt),t < (} under P is a (continuous) weak solution 
to (4.1) with initial distribution tt up to time Since aij,V G C^{D), weak uniqueness holds for 
solutions of (4.1) (see [33]). So under our conservativeness assumption, (4.1) has a conservative weak 
solution with initial distribution vr which is unique in distribution. By standard techniques (cf. the 
proof of [1, Proposition 1.2.1]), any weak solution to (4.1) gives rise to a solution to the martingale 
problem for {Q,T>(Q)). This implies that ( = oo because {X, K) with the initial distribution vr is a 
conservative solution to (4.1), by assumption. We conclude that tt is a stationary distribution for 
(4.1). □ 



We will present an easily verifiable condition on V for Theorem 4.1 to be applicable. The 
following preliminary result for diffusions without inert drift may have interest of its own. 

Recall that Wq'^{D) is the closure of the space C^{D) of smooth functions with compact 

/ ^ \l/2 

support in D under the Sobolev norm ||tt||i,2 := ( 

Theorem 4.2 Let D C M.'^ be a domain (i.e. a connected open set) and A{x) = {aij{x)) be a 
measurable dx d matrix-valued function on D that is uniformly elliptic and bounded. Suppose that 
if is a function in Wq''^{D) that is positive on D. Then the minimal diffusion X on D having 
infinitesimal generator C = 2^ j=i ^•^ conservative. 

Proof. The process X is the symmetric diffusion (with respect to the symmetrizing measure 
(/5(x)^ dx) associated with the Dirichet form {S,J^) in L'^ {D , ipipcf' dx), where 



£{u,v) = 
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and T is the closure of C^{D) with respect to the norm ^/Sl■, here we define 

£i{u,u) := £(u,u) + / u{x)'^ip{x)^ dx. 
Jd 

Let be the symmetric diffusion in D with respect to Lebesgue measure on D associated with 
the Dirichlet form {£^ ,Wq''^{D)) in L'^{D,dx), where 

I r d 

£^{u,v) = - / aij{x)diu{x)djv{x)dx. 

Then by [13, Theorems 2.6 and 2.8], X can be obtained from X^ through a Girsanov transform 
using the martingale dZt = ip{X^ )^^dMf , where M''' is the martingale part of the Fukushima 
decomposition for ip{X^ ) — ip{Xq). We now conclude from [13, Theorem 2.6(ii)] that X is conser- 
vative. □ 



Theorem 4.3 Suppose that T is a constant symmetric positive definite matrix. Let D C M be 
a bounded domain, a = o"(.x) be a d x d matrix that is on D so that A = g is uniformly 
elliptic and bounded on D. Suppose that V is a function in D such that (p = e^^/^ G Wq''^{D). 
Then the (minimal) solution of (4.1) with initial distribution tt is conservative and has tt as its 
stationary distribution. 

Proof. Let (X, P-j,) denote the symmetric diffusion process with infinitesimal generator 

1 

By Theorem 4.2, {X,¥-j;) is conservative. Let Kt = Kq - ^ J^{rV){Xs) ds. Define for t > 0, 

Mt := exp (^a-\Xs)Ks dBs-^J^ \a-\Xs)Ks\^ ds^ . 

Clearly {Mt,t > 0} is a continuous positive Pa;-local martingale for every x e D with respect to 
the minimal augmented filtration {J^t,t > 0} of X. For n > 1, let r„ := inf{t > : \Kt\ > 2"}. 

Since X is conservative, we have Too '■= limn^oo = oo, Px-a.s., for all x £ D. 

Let P be the distribution of the symmetric stationary diffusion X with initial probability dis- 
tribution cie^^ dx in D and K = 0. The existence of such a measure follows, for example, from 
Theorem 4.2 with F = 0. Let Kq have the Gaussian distribution C2 exp(— (F~^y, y)). Define a new 
probability measure Q on .Fqo by 

^ = Mt on J^t for every t>0. 
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By the Girsanov theorem, the process 

Wt -Bt- I a-\Xs)Ksds 
Jo 

is a Brownian motion under Q up to the time T^, so under Q, 
UXt = a{Xt) dWt + h{Xt) dt - ^{AVV){Xt) dt + Kt dt, 

\dKt = -^rvviXt)dt, 



for < i < Too. 



(4.7) 



Note that {X, K) has initial distribution vr. 

Conversely, given a solution {X, K) of (4.7) under the measure Q, the process X is a conservative 
>Cj;-diffusion under the measure 



exp (^-a-\X,)Ks dWs-\j^ \a-\Xs)Ks\'' d^ 



So to prove that the solution to (4.1) is conservative, it suffices to show that Q(Too = oo) = 1. 

We are going to show that {M(,t > 0} is in fact a positive P-martingale. This implies that 
Q(roo = oo) = 1 because q{T^ > t) = Ep [M^] = 1 for every t > 0. 

Let E := D and E'a = E \J {A} be the one-point compactification of E. Recall Q is 

defined by (4.2) with T>{Q) = C^{E). Since (4.4) holds, by the same argument as in the proof 
of Theorem 4.1, we deduce that tt is a stationary measure for some solution Q on ]D)(R+,^a) to 
the martingale problem for {Q,'D{Q)). Let {X,K) denote the coordinate maps on E^ and set 
C := ini{t > : {Xt,Kt) = A}. Then {{Xt,Kt),t < Q satisfies the SDE (4.7), and consequently, it 
has the same distribution as {{Xt,Kt),t < Too} under Q. 

Note that the matrix is bounded so there is a constant ci > such that 

|(7~"^(x)v| < ci\v\ for every x E D and v G M*^. 

Since under Q, Kt has the same Gaussian distribution for every t > 0, there exist C2 and cs such 
that if r < C2, then 



Ej 



exp r 



(^r\a-\Xs)Ks\^^ < E ^ [exp (rcf | | ^) ] < cg for every s > 0. 



By Jensen's inequality applied with the measure ^l[o,tol(*) with to G (0, C2] and the function 
we have 



E; 



exp ^a-\Xs)Ks 



ds 



E = 



< E; 



1 /"*° 



(to \a-''{Xs)Ks 



ds 



to 



I 

Jo 



exp to 



a-\Xs)k, 



ds 



I fto 
< — C2,ds = C3. 



*o Jo 
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Define Nt = J* a-\Xs)Ks clBs. This is a martingale with respect to P and its quadratic 
variation {N)t is equal to JjJ \a^^{Xs)Ks\'^ ds under both P and Q. Note that exp( — A''^ — ^{N)^) is 
a positive local martingale with respect to P and hence a P-supermartingale. Recall that is the 
lifetime for {X, K), i.e., {Xt, K) = A for i > Too, and that by convention, every function f on E is 
extended to a function on £'a by setting /(A) = 0. In particular, we have Nf = Nt^x^. We have 

g-ivt-i(iv>; 



-2Nt 



< 1. 



Using Cauchy-Schwartz, 



exp I / 1(7 ^{Xs)Ksf ds 



E« 



AN), 



<(EQ[e-^^*])'^'(EQ[eW*]) 



1/2 



< (^E( 

< I E; 



exp |c7 ^(Xs)Ksj^ 
exp^*|a-H^.)i?. 



1/2 
1/2 



As we observed in the previous paragraph, the last term is bounded \i t < C2. It follows from 
Novikov's criterion (see [31, Proposition VIII. 1.15]) that {Mj, t G [0, C2]} is a uniformly integrable P- 
martingalc. It follows that Q(Too > C2) = 1. Using the Markov property, we have Q(roo = 00) = 1. 
Consequently, tt is a stationary distribution for {X,K) of (4.1). This proves the theorem. □ 



The next corollary follows immediately from Theorem 4.3 and the fact that the solution of (4.1) 
depends in a continuous way in its initial starting point (xq, feo)- 

Corollary 4.4 Under the conditions of Theorem 4-3, for every {xo,ko) E D x W^, the minimal 
solution of (4.1) with initial value {xo,ko) is conservative. 

The following remark gives some sufficient conditions for if G Wq''^{D) and thus for the condition 
of Theorem 4.1 to hold with e~^ = ip"^. 

Remark 4.5 (i) Let W'^''^{D) := {/ G L^{D,dx) : V/ G L'^{D,dx)}. It is known (see, for 
example, [17]) that a function u G W^'^{D) is in Wq''^{D) if and only if its quasi-continuous 
version vanishes quasi-everywhere on dD. So in particular, if m G W^''^{D) and if 

lim u{x) = for every z G dD, 

a;— »2 

then u G Wo'^(-D). 
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(ii) For any positive bounded function u G W^'^{D) n if = e"^/" G W^''^{D) n C^{D) 
and so Tlicorcm 4.1 is applicable for y = 2/77, in view of Theorem 4.3. This follows because 
/ := exp(— l/(2x)) is bounded and Lipschitz continuous on [0,k] for every A: > 1. So by the 
normal contraction property (cf. [17]), V = f{u) G Wq^'^(D). In particular, this is the case if 
u is the first positive eigenfunction for the Dirichlet Laplacian on D (if such a ground state 
exists). 

(iii) The assumption that e~^/^ G C^{D) n Wq''^{D) is close to being optimal for ensuring that 
the solution to (4.1) never hits the boundary of the domain D. Consider the case n = 1, 
D = R+, (7 = 1, y : M+ ^ R+, and remove the inert drift so that dX = - \W{X) dt + dB. 
Define a function : R+ ^ R by 

^{x) = exp(F(y)) dy . 

Then one can check that ^{Xt) is a local martingale with quadratic variation greater or equal 
to 1. In particular, it has a positive probability of taking arbitrarily large values during any 
finite time interval. If V now diverges at sufficiently slowly so that exp(y) is still integrable, 
this easily implies that this diffusion has a positive probability of reaching in any finite time. 

Let us now consider the case where exp(F(x)) = C/x" near x = 0. In this case, ^e~^/^ is 
L^-locally integrable near if and only if a > 1, which is almost the range of parameters for 
which exp(y) is no longer integrable at 0, namely, a>\. 



The remainder of this section is devoted to the analysis of the case when both a (and hence A) 
and r are the identity matrix. We consider the SDE 

{dXt = - \VV{Xt) dt + Ktdt+ dBt , 
\dKt = -^VV{Xt)dt . 

It is possible to show under very weak additional assumptions that tt is the only invariant 
measure for (4.8). This is not obvious for example in the case where V{x) = f{d{x,dD)) for some 
function / : which has a singularity at and is such that /(r) =0 for r larger than some 

(small) constant e. On the set {x : d{x, dD) > e} x W^, the diffusion (4.8) has then a deterministic 
component dK = 0. In particular, this shows that (4.8) is not hypoelliptic and that its transition 
probabilities are not absolutely continuous with respect to Lebesgue measure. 

For any multi-index a = {ai, . . . , a^}, we define the vector field VaV by 

(V^Vix)),- ^'^'^("^ 

dXj^dX(y^ . . . dXfy^ 

Denoting by |a| the size of the multi-index, we furthermore assume that 
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Assumption 4.6 There exists £ D such that the collection {^aV{x*)}\a\>o spans all ofM!^. 

Remark 4.7 In the case where D is smooth and V is of the form V{x) = f{d{x,dD)) for some 
smooth function / diverging at 0, Assumption 4.6 is satisfied as soon as there exists a point on the 
boundary such that its curvature has full rank. 

We then have 

Proposition 4.8 // e'^/^ ^ C'^(^D) n Wo'^(-D) and Assumption 4.6 holds, then Tr{dx,dy) := 
cqId exp(— l^(a;) — \y\^)dxdy is the unique invariant measure for (4-8). 

Proof. Let us first introduce the following concept. Given a Markov operator V over a Polish 
space X, we say that V is strong Feller at x if V<p is continuous at x for every bounded measurable 
function <p: X ^ M.. The proof of Proposition 4.8 is then based on the following fact, a version of 
which can be found for example in [25, Thm 3.16], which is a consequence of well-known results 
from [14, Chapter 4] and [27, Section 6.1.2]: If {Vt}t>o is a Markov semigroup over X such that 

(i) There exists t > and x e X such that Vt is strong Feller at x, 

(ii) One has x G suppn for every invariant probability measure tt for the semigroup {Vt}, 

then the semigroup {Vt} can have at most one invariant probability measure. 

Denote now by Vt the Markov semigroup generated by solutions to (4.8). It is easy to check 
that Assumption 4.6 means precisely that the operator dt — G, where Q is the generator of (4.8), 
satisfies Hormander's condition [18, 28] in some neighborhood of {x^,,K,t) for every K e and 
every t > 0. Since for any (f) G ^^(R''), the map $(t, y, K) = [Vt4>) {y, K) is a solution (in the sense 
of distributions) of the equation (pt — ^) $ = 0, this implies that Vt<t> is C°° in a neighborhood of 
(x*, 0) for every t > 0. In particular, Vt is strong Feller at (x*, 0) for every t > 0. 

It now remains to show that the point (.x*, 0) belongs to the support of every invariant measure 
of (4.8). This will be checked by showing first that (ii) follows from the following two properties: 

(iii) for all y ^ X, there exists s{y) > such that x* G suppVs(y){y, ■ )■ 

(iv) for every neighborhood U of x^ there exists a neighborhood U' G U and a time Tjj > such 
that inf y^u' Vt{y,U) > ioi every t <Tu- 

The argument as to why (iii) and (iv) imply (ii) is the following: 

Fix an arbitrary neighborhood U of and an arbitrary invariant measure tt. By the definition 
of an invariant measure, one then has 

■k{U)= e-'Vtiy,U)7ridy)dt> / e-^'^y^+'^Vs(y)+tiy,U) dt7r{dy) 

= f f e-^<y^+M{y', U)Vsiy) (2/' dy') dt TTidy) > . 

Jx Jo Ju' 
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In the last inequality, we used the facts that Vs(^y^{y,U') > by (iii) and Vt{y',U) > by (iv). 
Since U was arbitrary, this inequality is precisely (ii). 

Since (iv) is satisfied for every SDE with locally Lipschitz coefficients, it remains to check that 
(iii) is also satisfied. For this, it is enough to check that, for every {x, K) G M^*^, there exists T > 
such that Vt{x, K;A)>0 for every neighborhood A of the point (x*, 0). In order to show this, we 
are going to apply the Stroock-Varadhan support theorem [33] , so we consider the control system 

x = -^VV{x) + K + u{t) , k = ~VVix), (4.9) 

where u: [0,T] — > R'' is a smooth control. The claim is proved if we can show that for every 
(xq, Kq), there exists T > such that, for every e > there exists a control such that the solution 
to (4.9) at time T is located in an e-neighborhood of the point (x*,0). Our proof is based on the 
fact that, since we assumed that V{x) grows to +oo as x approaches dD, there exists a collection 
of i points xi, . . . ,Xi {i > n) such that the positive cone generated by W{xi), . . . , W{xi) is all 
of M^. (This fact has already been noted in the paragraph following (4.3).) Fix now an initial 
condition (xq, -f^o)- From the previous argument, there exist positive constants ai, ... ^a^ such that 
ELi aiVF(xi) = -Kq. Fix now T = YlUi and consider a family : [0,T] of smooth 

trajectories such that: 

(i) there are time intervals Ii of lengths greater than or equal to aj — £ such that {t) = Xi for 
teli, 

(ii) ^^(0) = xq and X^iT) = x*, 

(iii) one has /[o,T]\U/i - ^• 

Such a family of trajectories can easily be constructed (take a single trajectory going through all 
the Xi and run through it at the appropriate speed). It now suffices to take as control 

u{t) = Mt) + \yv{x,{t)) + ^J^ VViXM) , 

and to note that the solution to (4.9) is given by (^X^{t),Ko — J^'W{Xs{s))) . This solution has 
the desired properties by construction. □ 

Remark 4.9 The paper [8] also shows uniqueness of the invariant measure for a mechanical system 
for which hypoellipticity fails to hold on an open set. However, the techniques used there seem to 
be less straightforward than the argument used here and cannot easily be ported to our setting. 

Remark 4.10 Through all of this section, the assumption that exp(— y/2) G Wq''^{D) could have 
been replaced by the assumptions that V £ C°°, limx^QD V {x) = +oo, J"^ exp(— F(x)) dx < oo, 
and there exists a constant C G M such that the relation 

AV{x) < |Vy(x)p + C(l + V{x)) (4.10) 
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holds for every x ^ D. This is because (4.10) ensures that V{x) + is a Lyapunov function for 
the diffusion (4.8) and so guarantees the non-explosion of solutions. 

5 Weak convergence 

In this section, we will prove that reflecting diffusions with inert drift can be approximated by 
diffusions with smooth inert drifts. 

Let -D be a bounded Lipschitz domain in M'^ and use Sjjix) to denote the Euclidean distance 
between x and D^. A continuous function x i— > 6(x) is called a regularized distance function to D'^ 
if there are constants C2 > ci > such that 

(i) ciSoix) < 6{x) < C25d{x) for every x € W^; 

(ii) S G C°°{D) and for any multi-index a = {ai, • • • , aa) with |q;| := Ylk=i ^k, there is a constant 
Ca such that 

9l°l(5(x) 



(9xi)"i • • • {dxd)°''' 



< CaSoix)^ for every x e D. 



The existence of such a regularized distance function S is given by [36, Lemma 2.1]. For n > 1, 
define 

, , I exp I 1 1 for X G -D, , , 

[+00 for X G D'^. 

Observe that Vn G C°°{D) and for every multi-index a, the a-derivative of e~^" at x G D tends to 
zero as X ^ dD. So e^^" G C°°(M'^). Observe also that as n ^ 00, e~^"(^) increases to loix). 

Suppose that A = a'^a is a uniformly elliptic and bounded matrix having entries aij on D, 
and p is a function on D that is bounded between two positive constants. 

We will consider the smooth potential approximations to the symmetric reflecting diffusion 
process on D associated with generator C = ■^V{pAV) but with an additional inert drift. More 
precisely, let F be a positive definite constant symmetric d x d-matrix. As before, we use n(x) to 
denote the unit inward normal vector at x G dD and define u(x) := A(x)n(x). We consider the 
following special case of the SDE (2.1): 

= a{Xt) dBt + h{Xt) dt + n{Xs) dLt + Kt dt , 
= Tn{Xt)dLt, 

where Xt €z D for every t > 0, and L is continuous and non-decreasing with Lt = Jq 1i)d{Xs) dLg- 
Here B \s a, d-dimensional Brownian motion and the drift b is given by b = (61 , • • • , 6^) with 

hk{x) := ^pYfi=idi{p{x)(Hk{x)). 

Let be the diffusion given by 

= (7(xi")) + b(xi")) - ^(yivy„)(xj"^) dt, 
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with initial probability distribution TTn{dx) := c„ exp(— V^(x))p(.T)l£)(a;) dx, where B^"^^ is a d- 
dimensional Brownian motion, and c„ > is a normalizing constant so that 7rn(K'^) = 1- Hence 
X^'^^ is the minimal diffusion in D with infinitesimal generator 

which is symmetric with respect to the measure TTn- Clearly e('°g''-^")/2 g C{p) n Wl^{p) and 
so by Theorem 4.2, is conservative and never reaches dD. In addition, is a symmetric 
diffusion with respect to the measure 7r„ in D and its Dirichlet form (f("),:r(")) in L2(D, 7r„(dar)) 
is given by (see [36, Lemma 3.5]) 

^^"^(/,/) ■■=\l aij{x)dif{x)djfix)7rn{dx) , 

:= {/ G L2(Z),7r„) : f (")(/,/) < oo} . 

Without loss of generality, we assTimc that X^"' is defined on the canonical path space Q := 
C{[0,(X)),M.'^) with Xf^\u!) = uj{t). On fi, for every t > 0, there is a time-reversal operator rt, 
defined for a; € by 

/c^(t-s), ifO<s<t, 
n(a;)(s) := <^ (5.3) 
I a;(0), 11 s >t. 

Let P„ denote the law of X^") with initial distribution 7r„ and let {.?^/"'' , i > 0} denote the minimal 
augmented filtration generated by X^^\ Then it follows from the reversibility of X^"-* that for every 
t > 0, the measure on ^^"^ is invariant under the time-reversal operator rt (cf. [17, Lemma 
5.7.1]). We know that 

^(n) _ ^(n) ^ ^(n) ^ ^ (b(XH) - i(^VK)(Xi"))) ds, t > 0, (5.4) 

where Adj:"^^ := a{xi"'^) dB^"^ is a square-integrable martingale. On the other hand, by the 
Lyons-Zheng's forward and backward martingale decomposition (see [17, Theorem 5.7.1]), we have 
for every T > 0, 

^(n) _ ^{n) ^ 1 ^^(n) _ 1 ^^^n) _ ^ W J „ foj. ^ g [q^ . (5.5) 

Since a is bounded, for every T > 0, the law of {M^^\t > 0} under P„ is tight in the space 
C([0,T],M'') and so is ^{M^'^ - M^\) o rr- It follows that the laws of {X^''\ M^''\ {M^^ - 
-^^r- ) ° '''t) under P„ are tight in the space C([0, T], M^'^) for every T > 0. On the other hand, by 
(5.4)-(5.5), 

^^^^^^^ '■= -\ /VVF.)(X(")) = -isf ) - \{bP - bP_,) oTT-f^ b(x(")) ds. 
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Thus the laws of (X("),M("), b(xi"^) ds, H^''\ S^")) under P„ are tight with respect to the 
space C([0, T], M"'^^) for every T > 0. By almost the same argument as that for [29, Theorems 
4.1 and 4.4], P„) converges weakly in C([0, T],M'^) to a stationary reflecting diffusion X in 

D with normalized initial distribution cp{x) dx on D. Passing to a subsequence, if necessary, we 
conclude that 

((XW,m("), /Qb(xi"^)(is, 5(")),F„) converges weakly in C([0, T], M'^'^) to a process 

((X,M, !^h{Xs)ds, H, B),F). 

Now we apply the Skorokhod representation (see Theorem 3.1.8 in [16]) to construct the pro- 
cesses (xH,M("), J^h{xi'^^)ds, SH) and {X,M, /ob(X,)(is, H, B) on the same proba- 
bility space (J^,.?=',P) in such a way that (X("),M("), /ob(xi"^)ds, S^) converges to 
(X, M, J"Qb(X5)(is, H, B) a.s., on the time interval [0,1] in the supremum norm. Therefore, in 
view of (5.4), 

Xt = Xo + Mt+ h{Xs)ds + H for every t > 0. 
Jo 

By the proof of [29, Theorem 6.1 and remark 6.2], if is a continuous process of locally finite 
variation. On the other hand, since X is reflecting diffusion in a Lipschitz domain, by [12] it admits 
a Skorokhod decomposition. That is, 

Xt = Xo+ [ a{Xs)dBt+ [ h{Xs)ds+ [ A{Xs)n{Xs)dLs, t > 0, 
Jq Jq Jo 

where B is a Brownian motion with respect to the minimal augmented filtration generated by Xf, 

n is the unit inward normal vector field on dD, which is well defined a.e. with respect to the surface 

measure a, and L is a positive continuous additive functional of X with Revuz measure a. By the 

uniqueness of Doob- Meyer decomposition, we have 

Ht= f A{Xs)n{Xs) dLs, fort>0. 
Jo 

Define for t >0, 

^t"^ ■= - /* VF„(xi")) ds and Kt := T f n(X,) dL,. 
^ Jo Jo 



Set 



and 



Zt = exp (^j'^ a-\Xs)Ks dB, - ^ \a-\Xs)K,\^ ds^ t> 0, 



= exp a-\x(-^)K(-^ dB^^ - \ ^ \o-\xf))Kf\^ d.) , t> 0. 

Let dQ = ZicflP and dQn = By the Girsanov theorem, under Q„, (XH,ifH) satisfies the 

following equation 

'dX^^^ = a{xi''^)dwt^ + b(xf )) - l{AVVn)ixi''^)dt + i^f ^ dt, 

irvvnixi' 
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where W^""^ is a d-dimensional Brownian motion. On the other hand, by the proof of Theorem 2.1, 
{X, K) under Q is Brownian motion with inert drift, and satisfies 

[dXt = a{Xt) dWt + HXt) dt + (An)(Xt) dLt + Kt dt, 
\dKt = Tri{Xt)dLt, 

where is a d-dimensional Brownian motion. 

In the following, the initial distribution of {x'^\k^^) and {Xq,Kq) under P are taken to be 

c„lD(a;)/9(x)e-^"(^)-(^"'2/'2/) dxdy (5.6) 

and cl/j(a;)p(x)e~^^ ^'^''^^dxdy, respectively, where > and c > are normalizing constants 
chosen to make sure that the corresponding measures are probability measures. Since e('°sp-^n)/2 ^ 
C (D) DWq''^ (D) , we know from Theorem 4.3 that CnlD{x)p{x)e~^"^^^~^^ ^2/,2/) dxdy is a stationary 
measure for {X^'"\K^^^) under Q„. 

Theorem 5.1 For every T > 0, the law of under Q„ converges weakly in the space 

C([0,T],M2<i) to that of the law of{X,K) under <Q. 



Proof. Without loss of generality, we take T = 1. 



Observe that ^(X^) converges to A{X) uniformly on [0, 1] and Kf' = T A~^{Xr')dHr' ■ 
Let \H^^^^\f and ji^^^^jf denote the total variation process of and over the interval [0, t] , 

respectively. Then there is a constant c > independent of n > 1 such that 

It < It for every t > 0. 

On the other hand, we know from the proof of Theorem 6.1 in [29, p. 58-59] that 



sup E I 

n>l 



< oo for every t > 0. 



Hence by Theorem 2.2 of [22], converges to K in probability with respect to the uniform 

topology on C[0, 1]. This yields, by [22, Theorem 2.2] again, that /J cT-i(xi"V<^"^ f^^i""* ^ 
j} a^^{Xs)KsdBs as n ^ oo, in probability. By passing to a subsequence, we may assume that 

(n) 

the convergence is P-almost sure. We conclude that lim„^oo = Zi, P-a.s. 

Let $ be a continuous function on (C[0, 1])^ with < ^> < 1. Since ^ <i>iX,K), 

(n) 

P-a.s., and Z^ ' ^ Zi, P-a.s., by Fatou's lemma, 

Ev[^{X,K)Zi] < liminf Ep[$(x("),i^("))zJ"^] < limsupEp[$(x("\ i^("))zj"^] (5.7) 



and 

Ep[(l - ^){X,K)Zi] < liminfEp[(l - $)(x("), k("))z{"^]. (5.8) 
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Summing (5.7) and (5.8) we obtain Ep[Zi] < lim sup^^Q^^ E p [Zj" ] . Note that by the proof of 
Theorem 4.3, Z*^") is a continuous non-negative P-martingale, while by the proof of Theorem 2.1, 
M is a continuous P-martingale. Hence, Ep[Zi] = 1 = Eplzj""*] and, therefore the inequalities in 
(5.7) and (5.8) are in fact equalities. It follows that 

lim Q„[$(X('^),i^("))] = lim Ep[$(X('^), ii:("))z|")] = Ep[$(X, K)Zi] =Q[$(X,K)]. 

n— ♦oo n-^oo 

This proves the weak convergence of {X^'^\ K^"^^) under Q„ to {X,K) under Q. □ 

Theorem 5.2 Consider the SDE (5.2) on a hounded Lipschitz domain D C M"^. A stationary 
distribution for the solution {X, K) to (5.2) is 

Tr{dx,dy) = cilD{x)p{x)e~^^ ^'^^ dxdy, 

where ci is the normalizing constant so that it{D x R'') = 1. 

Proof. Recall the notation from previous proofs in this section. We start {X^'^\K^'^^) with the 
stationary distribution 7r„ of Theorem 4.3, namely, 

Q„((xi"\ i^i")) eA)= 7r,(yl) = / e-^"(^V(x)e-(r"'s''^) dxdy, AcDx R'^, 

J A 

where c„ is a normalizing constant. Note that here we have e~^"^^^p{x) in place of e"^^^^ in 
Theorem 4.1. Clearly 7r„ converge weakly on D x M*^ to the probability tt. Also, the laws of 
converge weakly to the Q law of {X,B,K). If / is any continuous and bounded 
function on R^'', then for to < 1, 

EQf{Xt„Kt,)= lim EQ„/(x(;),i^(;))= lim / f{x,y)n,,idx,dy) 
= /_ f{x,y)Tr{dx,dy). 
This shows that the tt is a stationary distribution for the solution to (2.1). □ 

6 Irreducibility 

In this section, D is a bounded domain in R'^. We will prove uniqueness of the stationary 
distribution for (X, K) under assumptions stronger than those in previous sections. Let Qx,y denote 
the distribution of {X, K) with inert drift starting from (x, y). We say that {X, K) is irreducible in 
the sense of Harris if there exists a positive measure ji on D xW^ and to > such that if //(A) > 0, 
then for all (x, y) eD x R'*, Q^,j,((Xto, i^tj e A) > 0. 

Let r be a positive definite constant symmetric d x d-matrix. 
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Theorem 6.1 Assume that a is the identity matrix (consequently u = nj, b = i^^logp, and p is 
a function on D that is bounded between two positive constants. Then the solution {X, K) to 
(5.2) is irreducible in the sense of Harris. 

Proof. We first assume p = 1 on D. Let X denote a solution to (2.2), i.e., the usual normally 
reflecting Brownian motion in D, let L be the local time of X on dD, and Kf = yo + n{Xs) dLg. 
The distribution of (X, K) starting from (Xq, Kq) = (x, y) will be denoted P^,?/- First, we will prove 
Harris irreducibility for (X, K) under Pa;,j/, with respect to 2d-dimensional Lebesgue measure. 

In order to do so, our main ingredient will be that for i > the law of {Xt, Kt) has a component 
that has a strictly positive density with respect to Lebesgue measure on some open set. The idea 
behind the proof of this result is that one can find d small balls {Bj}j^-^ on the boundary dD such 
that 

(i) The normal vector n(a;) is 'almost' parallel to the j-th unit vector for x e Bj. 

(ii) With positive probability, the process Xg has visited all of the Bj^s in chronological order 
before time t, but has not hit the boundary dD anywhere outside of the i?/s. 

(iii) The amounts Sj of local time that X has accumulated on the Bj's are 'almost' independent. 

This suggests that at time t, the law of Kt has a positive component which 'almost' looks like the 
law of a random vector with independent components, each of them having a density with respect 
to Lebesgue measure. It follows that the law of Kt has a component which has a density with 
respect to d-dimensional Lebesgue measure. Since, as long as X is in the interior of the domain D, 
it is just Brownian motion with drift and K remains constant, we conclude that the law of (Xt, Kt) 
has a density with respect to 2£i-dimensional Lebesgue measure. 

The detailed proof is broken into several distinct steps. In the first step, we use a support 
theorem and a controllability argument to show that given any point zi E D, e > and any final 
time to, the process {Xt,Kt) has a positive probability to be in an e-neighborhood of {zi,0) after 
time to/2, whatever its initial condition. In the second step, we present a review of excursion theory 
and show how the path of a reflecting Brownian motion can be decomposed into a collection of 
excursions and how reflecting Brownian motion up to the first hitting time of some subset U C D 
can be constructed from the excursions of a reflecting Brownian motion conditioned never to hit 
U by adding a 'last excursion' after a suitably chosen amount of local time spent at the boundary. 
This construction is then used in the third step to 'stitch together' a reflecting Brownian motion 
from d independent reflecting Brownian motions . In the final step, we show how to condition 
each of the Y^'s on hitting the boundary dD only in Bj and deduce from this that the local time 
Kt has a density with respect to Lebesgue measure. We conclude by showing how to combine these 
results to obtain the desired Harris irreducibility. 

Step 1. Fix any to,r > and zi G D. In this step, we will show that for any {xo,yo) E D x 
there exists pi > such that P^o,j,o((Xt^/2, -^to/2) € B{zi,r) x B{0,r)) > pi. 
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We recall the deterministic Skorokhod problem in D with normal vector of reflection. Suppose 
a continuous function / : [0, T] is such that /(O) € D. Then the Skorokhod problem is to find 

a continuous function g : [Q,T] D and a non-decreasing function I : [0,T] [0, oo), such that 
^(0) = 0, g{0) = /(O), /(f d£, = 0, and g{t) = f{t) + n{g{s)) Ms- It has been proved in 

[24] that the Skorokhod problem has a unique solution {g, t) in every domain. 

Since D is a bounded smooth domain, the set {n(a:)/|n(x)|, x G dD} is the whole unit sphere 
in W^. Find x\ G dD such that n{x\) = — coyo for some cq > 0. It is elementary to construct a 
continuous function / : [0,to/2] ^ M'^ such that /(O) = xq, f{t) e D ior t e (0,to/4), /(to/4) = xi, 
f is linear on [to/4, 3to/8], f{Mo/8) = xi + yo, f{t)-yo G for t G (3to/8, to/2), and /(to/2)-yo = 
Zi. It is straightforward to check that the pair (5,^) that solves the Skorokhod problem for / has 
the following properties: g{t) = f{t) for t G (0,to/4), g{t) = xi for t G [to/4, 3io/8], g{t) = /(t) -yo 
for t G (3to/8, to/2), and J^"^^ n{g{s)) d4 = -yo- _ 

For a function : [0,to/2] ^ M.^ with /^(O) G D, let (5\^^) denote the solution of the 
Skorokhod problem for f^. Let Bc{f,S) be the ball in C([0, to/2], R*^) centered at /, with radius 
6, in the supremum norm. By Theorem 2.1 and Remark 2.1 of [24], for any r > there exists 
6 G (0,r/2), such that if G Bc{f,S), then c/^ G Bc{g,r/2). This implies that 

fto/2 ^ ^ /•to/2 

sup ' ' 



te[o,to/2] 



/•to// /•to// 
/ n{g\s))Ml- n{g{s))di, 
Jo Jo 



< sup {\f\t)-m\ + \g\t)-gm 

te[o,to/2] 

< ^ + r/2 < r. 

Thus, if /I G Bc(/,5), then (5i(to/2), Jq"/' n((7i(s)) d^j) G B{zi,r) x ^(-2/o,r). Let denote 
the distribution of standard Brownian motion. By the support theorem for Brownian motion, 
ho{Sc{f,S)) > Pi, for some pi > 0. Hence, P^o,2/o((^to/2>^to/2) e B{zi,r) x B{0,r)) > pi > 0. 

Step 2. This step is mostly a review of the excursion theory needed in the rest of the argument. 
See, e.g., [26] for the foundations of excursion theory in abstract settings and [9] for the special 
case of excursions of Brownian motion. See also [19] for excursions of reflecting Brownian motion 
on domains. Although [9] does not discuss reflecting Brownian motion, all the results we need 
from that book readily apply in the present context. We will use two different but closely related 
exit systems. The first one represents excursions of refiecting Brownian motion from dD. 

We consider X under a probability measure P^;, i.e., X denotes reflecting Brownian motion 
without inert drift. 

An exit system for excursions of reflecting Brownian motion X from dD is a pair (L^,Ha.) 
consisting of a positive continuous additive functional and a family of excursion laws {H.x}xedD- 
We will soon show that = Lf. Let A denote a "cemetery" point outside and let C be the space 
of all functions / : [0, 00) — U {A} which are continuous and take values in R*^ on some interval 
[0, C), and are equal to A on [C,oo). For x G dD, the excursion law Hx is a cr-finite (positive) 



32 



measure on C, such that the canonical process is strong Markov on (to,oo) for every to > 0, with 
the transition probabihties of Brownian motion killed upon hitting dD. Moreover, H^- gives zero 
mass to paths which do not start from x. We will be concerned only with "standard" excursion 
laws; see Definition 3.2 of [9]. For every x G dD there exists a standard excursion law Hx in D, 
unique up to a multiplicative constant. 

Excursions of X from dD will be denoted e or eg, i.e., if s < u, Xs,Xu G dD, and Xt ^ dD for 
t G {s,u), then = {es{t) = Xf^g, t ^ [^jU — s)} and Ci^s) = u — s. By convention, es{t) = A for 
t>C- So et = A if inf{r > t : G dD} = t. Let £u = {cs : s < u}. 

Let at = inf{s > : L* > t} and let / be the set of left endpoints of all connected open 
components of (0, oo) \ {t > : G dD}. The following is a special case of the exit system 
formula of [26]. For every x £ D, 



.tei 



E, 



/•OO "I r /"OO 

/ Z,^Uxi^^){f)ds =Ex / ZtUxMldL; , (6.1 
Jo J L-'o 



where Zt is a predictable process and f : C ^ [0, oo) is a universally measurable function which 
vanishes on those excursions et identically equal to A. Here and elsewhere H^(/) = fdHx- 

The normalization of the exit system is somewhat arbitrary, for example, if {Lf,Hx) is an exit 
system and c G (0, oo) is a constant then {cL^ , {1 / c)Ux) is also an exit system. Let Fy denote 
the distribution of Brownian motion starting from y and killed upon exiting D. Theorem 7.2 of 
[9] shows how to choose a "canonical" exit system; that theorem is stated for the usual planar 
Brownian motion but it is easy to check that both the statement and the proof apply to reflecting 
Brownian motion in D C M*^. According to that result, we can take to be the continuous additive 
functional whose Revuz measure is a constant multiple of the surface area measure dx on dD and 
the Hj;'s to be standard excursion laws normalized so that for some constant ci G (0, oo). 



H^(^) 



1, 



cilim^Pf^,„(.)(A) 



(6.2) 



for any event A in a cr-field generated by the process on an interval [to, oo) for any to > 0. The 
Revuz measure of L is C2 dx on dD. We choose ci so that (Lj, H^:) is an exit system. 

We will now discuss another exit system, for a different process X'. Let U C D he a fixed 
closed ball with positive radius, and let X' be the process X conditioned by the event {T^ > ai}, 
where = M{t > : X G [/}. One can show using Theorem 2.1 and Remark 2.1 of [24] that for 
any starting point in D\U, the probability of {T^ > is greater than 0. It is easy to see that 
{X[,Lt) is a time-homogeneous Markov process. For notational consistency, we will write (X^,LJ) 
instead of {X[,Li). 

We will now describe an exit system (L^, ^) for {X[, L[) from the closed set D x [0, oo). We 
will construct this exit system on the basis of {Lt, Hx) because of the way that X' has been defined 
in relation to X. It is clear that L' does not change within any excursion interval of X' away from 
dD, so we will assume that ^ is a measure on paths representing X' only. The local time L' is 
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the continuous additive functional with Rcvuz measure C2 dx on dD. For ^ > 1 wc let ^ = H^;. 
Let ¥y denote the distribution of Brownian motion starting from y ^ D\U, conditioned to hit dD 
before hitting U, and killed upon exiting D. For i < 1, we have 

H;,(A)=cilimlpf^,„(,)(^). (6.3) 

Let ^* C C be the event that the path hits U. It follows from (6.2) and (6.3) that for £ < 1, 

H^,,(^) = H,(^\A,). (6.4) 

One can deduce easily from (6.2) and standard estimates for Brownian motion that for some 
C3, C4 G (0, oo) and all x G dD, 

C3 < H^(A*) < C4. (6.5) 

Let a'-i- = inf{s > : > t}. The exit system formula (6.1) and (6.4) imply that we can construct 
X (on a random interval, to be specified below) using X' as a building block, in the following 
way. Suppose that X' is given. We enlarge the probability space, if necessary, and construct a 
Poisson point process S with state space [0, oo) x C whose intensity measure conditional on the 
whole trajectory {X^, t > 0} is given by 

/•1AS2 

/x([si,S2] X F) = / Ilx',{FnA^)dt. (6.6) 

Since /i([0, oo) x C) < oo, the Poisson point process E may be empty; that is, if the Poisson process 
is viewed as a random measure, then the support of that measure may be empty. Consider the case 
when it is not empty and let Si be the minimum of the first coordinates of points in £. Note that 
there can be only one point (5i, e^J G £ with first coordinate S\, because of (6.5). By convention, 
let = oo if = 0. Recall that = inf{t > : Xt e U} and let 

Tf = mf{t >0:X^eU}, 
% = a's^ + inf{i > : esi {t) G U}. 

It follows from the exit system formula (6.1) that the distribution of the process 



'x[ ifo<t<rf a4^. 



^e5i(i-(^5j iff /0anda^^ <i<r*, 
is the same as the distribution of {Xt, < t < T^}- 

Step 3. We will now construct reflecting Brownian motion in D from several trajectories, including 
a family of independent paths. 

Let Uj = B{zj, r) for j = 1, . . . , d + 1, where zj G D and r > are chosen so that Uj HUk = $ 
for j ^ k, and Ui<j<d+i Uj C D. 
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Recall from the last step how the process X was constructed from a process X' . Fix some 
xi G Ui and let X^ be a process starting from Xq = .xi, with the same transition probabilities as 
X' , relative to U2- We then construct 1"^ based on X^, by adding an excursion that hits U2-, in 
the same way as X was constructed from X' . We thus obtain a process {1^^,0 < t < Ti}, where 
Ti = 'mi{t > Q : Y^ ^ U2} 1 whose distribution is that of reflecting Brownian motion in D starting 
with the uniform distribution on Uj, observed until the first hit of ?72- 

We next construct a family of independent reflecting Brownian motions {Y^}i<j<d- For a 
fixed j = 2, . . . ,d, we let X^ be a process with the same transition probabilities as X', relative to 
Uj-^-i, and initial distribution uniform in Uj. We then construct based on X^ , by adding an 
excursion that hits t^j+i, in the same way as X was constructed from X' . We thus obtain a process 
{y/,0 < t < Tj}, where Tj = inf{t > : Y^^ £ Uj+i}, whose distribution is that of reflecting 
Brownian motion in D, observed until the first hit of t/j+i. 

Note that for some C5 > and all x,y e Uj+i, j = I, . . . ,d, 

^x{Xi e dy and Xt ^ dD for t G [0, 1]) > cg dy. 

We can assume that all X-^'s and Y-^ 's are defined on the same probability space. The last formula 
and standard coupling techniques show that on an enlarged probability space, there exist reflecting 
Brownian motions , j = 1, . . . ,d, with the following properties. For l<j<d — 1, Zq = Yj,^, 
and for some cg > 0, 

P (^Z( = and Zi i dD for t G [0, 1] | {Y%<k<j,{Z%<k<,_i^ > cq. (6.7) 

The process Z^ does not depend otherwise on {i^'^}i<jfc<d and {Z^}k^j. We define Z'^ as a reflecting 
Brownian motion in D with Zq = Yj.^ but otherwise independent of {l^'^}i<jfc<d and {.^'^}i<jfc<d-i- 
Let 

Fj = \z{ = Y^-^^ and ^ dD for i G [0, 1]} . 

We define a process X* as follows. We let X^ = Y^ io\ Q < t < T\. If F{ holds, then we let 

XI = Z}^j.^ for t>Ti. If Fi holds, then we let X^ = Z}_j,_^ for te[Ti,Ti + 1] and X^ = YI^j,^_^ 
for t G [Ti + l,Ti + 1 + T2]. We proceed by induction. Suppose that X^ has been defined so far 
only for 

t G [0,Ti + l + r2 + l + ---+Tfc], 
for some k < d. If holds, then we let 

— ^t-Ti-l-T-z-l Tfe 

for t > Ti + 1 + T2 + 1 + • • • + Tfc. If Ffe holds, then we let 

— ^t-Ti-i-r2-i Tfc 
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for t G [Ti + 1 + T2 + 1 + • • • + Tfe, Ti + 1 + Ta + 1 + • • • + Tfe + 1] and 

y* vk+1 

— ^ t-Ti-l-T2-l 

for t G [Ti + 1 + T2 + 1 + • • • + Tfe + 1, Ti + 1 + T2 + 1 + • • • + Tfc + 1 + Tk+i]. We let 

— ^t-ri-i-T2-i Td 

for t > Ti + 1 + r2 + 1 + • • • + Trf. 

By construction, X* is a reflecting Brownian motion in D starting from xi. Note that in view 
of (6.7), conditional on > 0}, j = l,...,d, there is at least probability Cg that X* is a 

time-shifted path of on an appropriate interval, for all j = 1, ... , d. 

Step 4- In this step, we will show that with a positive probability, the process K can have "almost" 
independent and "almost" perpendicular increments over disjoint time intervals. Moreover, the 
distributions of the increments have densities in an appropriate sense. 

We find d points yi,... ,yd G dD such that the n(?/j)'s point in d orthogonal directions. Let 
Cj = {z eW^ : Z{n{yj), z) < So}, for some > so small that for any zj G Cj, j = 1, . . . ,d, the 
vectors {zj} are linearly independent. Let (5i > be so small that for every j = 1, . . . ,d, and any 
y G dD n B{yj,5i), we have n(y) G Cj. 

Let be the local time of X^ on dD and = inf{s > : > t}. It is easy to sec that 
there exists P2 > such that with probability greater than p2, for every j = 1,. . . ,d, we have 
xi ^dD\ B{yj, 5i), for t G [0, (t{]. Let 

Rj = sup{i < Tj : y/ G dD} and Sj = L^j^. . 

Let be the event that for every j = 1, . . . ,d, Xj ^ dD \ B{yj,6i) for t G [0, a{] and Sj < (t-[. 
Then (6.5) shows that ^x^{F*) > P2{1 - e'"'')'^. 

Let Kf = r /* n(X"'j ) ds and note that if holds, then G TCj for all j = 1,. . . ,d and 
t G [0, 1]. Define for any < < 6^ < A; for A; = 1, • • • ,d, 

A([ai, [a2, 62], . . . , [a^, 6^]) = {K^^ +KI + --- + : tk G [a^, for 1 < A; < d}. 

It is easy to show using the definition of C/s that the d-dimensional volume of A([ai, 61], . . . , [a^, ba]) 
is bounded below by C7 ]^^<j;,<^(6jfc — a^), and bounded above by C8ni<jt<d(^ik ~ o,k)- 
Let us consider the processes defined above, conditioned on the a-field 

g = a(^{Ki,te[0, l]}i<,<d). 

By (6.5) and (6.6), conditional on Q, the random variables Sj = L-'j^^ , j = 1, • ■ • ,d, have distributions 
whose densities on [0, 1] are bounded below. In view of our remarks on the volume of A, it follows 
that conditional on Q, the vector Kl + • • • + Kg has a density with respect to cZ-dimensional 
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Lebesgue measure that is bounded below by cg > on a ball f7* with positive radius. We now 
remove the conditioning to conclude that Kl.^ + • • • + K'^.^ has a component with a density with 
respect to d-dimcnsional Lebesgue measure that is bounded below on U^. 



Define = T Jq n(X*s)(iL<j and T* = J2'j=i'^j^ where L is the boundary local time for 
reflecting Brownian motion X*. Using conditioning on F^, we see that the distribution of K^^ 
has a component with density greater than cg on i7*. Since K* does not change when X* is 
inside the domain and X* is a reflecting Brownian motion, we conclude that {X^^_^-^, K^^^^^) has a 
component with a density with respect to 2(i-dimcnsional Lebesgue measure on a non-empty open 
set. It follows that for some fixed t^, > 0, (X^* ) has a component with a density with respect 
to 2(i-dimensional Lebesgue measure on a non-empty open set. 

We leave it to the reader to check that the argument can be easily modified so that we can 
show that for any fixed to > 0, {X*^^^^ ^ component with a strictly positive density with 

respect to 2d-dimensional Lebesgue measure on a non-empty open set. We can now combine this 
with the result of Step 1 using the Markov property to see that for some non-empty open set U and 
any starting point {Xo,Kq) = {xo,yo), the process {XtQjKtg) has a positive density with respect to 
2(i-dimensional Lebesgue measure on U under ^xo,yo- 

By the Girsanov theorem the same conclusion holds for (X, K) under the measure Qxo,yoj which, 
by the proof of Theorem 2.1, is the reflecting diffusion with inert drift. This proves the theorem in 
the case p = 1. 

Now let us consider the general case when p G C'^{D) is bounded between two positive constants. 
In this case, {X,K) is a solution to the following special case of (5.2): 

f dXt = dBt + \V log p{Xt) dt + n{Xt) dt + Kt dt, 
\dKt = rn{Xt)dLt. 

Let {{X, K), P} denote the solution to (6.8), and let {(X°, K^), P°} correspond to the special case 
/9 = 1, with the same initial distribution. Note that our assumptions on p imply that Vlogp G 
C^{D). Hence P and P*^ are related by Girsanov transform 



^=exp(^- V\ogp{X^)dBs--J^ \V log p{X^)\'ds 



on J='t. 



The Harris irreducibility of {X,K) now follows from that for (X ,K ). □ 



Theorem 6.2 Consider the SDE (5.2) with a being the identity matrix, b = 2 log p, u = n and 
p € C'^{D) being bounded between two positive constants. The probability distribution 7:{dx,dy) 
defined by 

7r{A)= I cilD{x)p{x)e-^^~^y^y^ dxdy, A(zDxW^, 

J A 

is the only invariant measure for the solution {X,K) to (5.2). 
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Proof. In view of Theorem 5.2, all we have to show is uniqueness. If there were more than one 
invariant measure, at least two of them (say, /i and i/) would be mutually singular by Birkhoff 's 
ergodic theorem [32]. However, we have just shown that there exists a strictly positive measure 
ip which is absolutely continuous with respect to any transition probability, so that in particular, 
ip <^ jjL and il) <^ V. Since /x _L by assumption, there exists a set A such that n{A) = and 
v{A^) = 0. Therefore, one must have il^{A) = il){A'^) = which contradicts the fact that the 
measure is non-zero. □ 
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